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Motivated by the physics of graphene, we consider a model of N species of 2+1 
dimensional four-component massless Dirac fermions interacting through a three di- 
mensional instantaneous Coulomb interaction. We show that in the limit of infinitely 
strong Coulomb interaction, the system approaches a quantum critical point, at 
least for sufficiently large fermion degeneracy. In this regime, the system exhibits 
invariance under scale transformations in which time and space scale by different 
factors. The elementary excitations are fermions with dispersion relation cu ~ p z , 
where the dynamic critical exponent z depends on N. In the limit of large N, we find 
z = 1—4/ ( / k 2 N)+0(N^ 2 ). We argue that due to the numerically large Coulomb cou- 
pling, graphene (freely suspended) in vacuum stays near the scale-invariant regime 
in a large momentum window, before eventually flowing to the trivial fixed point at 
very low momentum scales. 

PACS numbers: 73.63.Bd, 05.10.Cc 

I. INTRODUCTION. 



The physics of graphene, which has recently been realized in experiment has currently 
attracted considerable interest. A graphene sheet is a two dimensional (2D) hexagonal lattice 
of carbon atoms, and its defining feature is the existence of two special points in the Brillouin 
zone around which the electron energy, according to band-structure calculations, should have 
a linear dependence on its momentum ("Dirac cones"), as in relativist* theories Q. Many 
interesting behaviors, for example, the quantum Hall effect at half-integer filling fractions, 
are attributed to the quasirelativistic behavior of the low energy excitations [1, Q, Ej] . 

Graphene, on the other hand, differs from relativistically invariant systems in one crucial 
aspect. Namely, the electromagnetic interaction between fermion quasiparticles is mediated 
by the photon whose velocity c is practically infinite compared to the fermion velocity 
v ~ c/300. The interaction therefore is an instantaneous Coulomb repulsion which breaks 
relativistic invariance. 

As we shall see later, the importance of the Coulomb interaction is controlled by the 
parameter 

e 2 N 
16e hv ' 

where N = 2 is the spin degeneracy. It is similar to the fine structure constant, but the 
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speed of light has been replaced by the Fermi velocity v. As the result, for graphene in 
vacuum this parameter is around 3 or 4 (it is reduced if graphene resides on a substrate with 
a large dielectric constant). The largeness of A makes unreliable any calculation based on 
the simple perturbation theory in the interaction strength, and seems to indicate that real 
graphene is hopelessly beyond quantitative theoretical control. 

In this paper, we show that in the idealized limit A — > oo the system becomes, in a certain 
sense, simple again. We shall argue that, at least for sufficiently large N, and perhaps even 
for iV = 2,A^ooisa limit where the system is tuned to quantum criticality. 

Quantum critical points (sj play an important role in condensed matter physics. In 
many cases, they are described by relativistically invariant conformal field theories. The 
relativistic invariance of these theories is not connected to the relativity of space-time, but 
is a manifestation of an emergent Lorentz symmetry. In our case, the infinitely strong 
Coulomb interaction destroys the relativistic invariance. Instead, we shall see that the 
quantum critical point A = oo is characterized by a nontrivial dynamic critical exponent z, 
whose value is computable in the large- N limit, 

z = 1-^ + 0^). (2) 

The elementary excitations are fermions with dispersion relation to ~ p z , instead of the linear 
dispersion uj = vp. The Dirac cones are therefore replaced by "Dirac cusps." Moreover, the 
full dynamics is invariant under the scale transformation 

t 1% x^lx (3) 

instead of the usual "relativistic" scale transformation 

t — > It, x — > Ix. (4) 

The limit A — > oo can therefore be considered as a nice idealization of graphene. While 
it is only an idealization, we think that its simplicity justifies our study. Moreover, one 
may hope that once the limit A — ► oo is understood, the case of finite but large A can be 
accommodated by treating A" 1 as a small coefficient of a relevant deformation that takes 
the system away from the critical point. 

The structure of this paper is as follows. In Sec. [Ill we describe our model. We write down 
the renormalization group (RG) equation in Sec. DTD and discuss the strong-coupling fixed 
point where the scaling behavior ([3]) is realized. The running of the fermion velocity at finite 
Coulomb coupling is considered in Sec. [TV] In Sec. |VJ we discuss finite iV and the relevance 
of our model to real graphene. We conclude in Sec. IVI1 While a large fraction of technical 
calculations presented in this paper are not new, we believe that the strong-coupling limit 
has not been previously identified as a quantum critical point. 
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II. MODEL 

We shall consider a system of N (2+1) dimensional [(2+l)D] four- component massless 
Dirac fermions with velocity v interacting through an instantaneous three dimensional (3D) 
Coulomb interaction. The Euclidean action is (in this paper, we set h = 1) 

N _ _ _ . 

S = -J2 Utd 2 X (iJ a J°d ^a + Vij a fd^ a + iA ^ a J ^ a ) + — dt d*X (^A ) 2 . (5) 

0=1 ^ y •* 



Our notations are as follows. The fields ip a are four-component fermion fields, and a labels 
different species of fermions. In real graphene N, is equal to 2 due to the spin degeneracy. 
The 7's are (2+l)D Dirac matrices satisfying {7^,7"} = 25 tJ,v , and can be chosen as 4 x 4 
matrices, e.g., 7 = a 3 ® a 3 , 7* = a 1 <S> 1. Each fermion can be thought of as a pair of 
two-component fermions with opposite parities 0, 0]- A§ is the Coulomb potential. The 
action (jSJ) contains a (2+1 )D part, which contains the kinetic term for the fermion and the 
interaction between the fermion and the Coulomb potential, and a (3+1) dimensional part, 
which is the kinetic term for the Coulomb potential. A more general model was considered 
in Ref. 0] in the context of the quantum Hall effect. If graphene is in vacuum (as it is the 
case for freely suspended graphene sheets [l(|), then g 2 = e 2 /e , where e is the electron 
charge and 6q is the vacuum permeability. In the presence of a substrate with a dielectric 
constant e, the effective charge is reduced, 



e 2 



2 = _^^ 

l + ee 1 ; 

In Eq. (JSJ), we have not included any contact four- Fermi interactions, which are irrelevant 
at weak coupling. At strong coupling, however, these interactions develop nontrivial fixed 



points [llj . We shall assume that these four- Fermi interactions start out with small couplings 
and flow to the trivial fixed point. 

We will be particularly interested in the limit of infinite Coulomb repulsion g 2 — > 00. In 
this limit, there is no kinetic term for Aq in the bare action (jSJ); however, an effective kinetic 
term will be generated by the fermion loop. We find it useful to keep g 2 large but finite for 
the purpose of regularization and for the discussion of the real graphene. 

Without the coupling to the scalar potential A , the theory is that of free Dirac fermions 
which is invariant under the relativistic scale transformation (141) . Our goal is to see that 
after coupling to the scalar potential the system remains scale invariant, but under a more 
general type of scale transformation with a (generally) fractional z. 



A relativistic counterpart of our model was considered previously. In Refs. [12], [13j free 
fermions are coupled to a gauge field A^. At sufficiently large N, the infrared limit of the 
new system is described by a conformal field theory. (For iV smaller than some critical 
value, which is still not exactly known, it was argued that the system develops a mass 



gap [14j.) The difference between our case and the case considered in Refs. 12|, |l3( is the 



lack of Lorentz invariance due to the instantaneous Coulomb interaction. The fact that in 
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graphene Coulomb interaction induces a logarithmic renormalization of the fermion velocity 



15 . To make the 



and leads to logarithmic corrections in thermodynamics is well known 
paper self-contained, we repeat some of the calculations in the literature. To start, let us 
state the Feynman rules that follow from Eq. (jHJ), 

(i) The fermion propagator is 

GoiP) = i = K- (7) 

Here we use quasirelativistic notation, where p stays for (po,p) (for example d 3 p = 
dp dp) and = 7V0 + v l 4 Pi P 2 = Vo + w 2 |p| 2 , and p is the 2D momentum vector. 

(ii) The Aq propagator is the integral of the 3D propagator g 2 /(p 2 z + \p\ 2 ) over the mo- 
mentum component perpendicular to the 2D plane, p z , 

and is simply the 2D Fourier transform of the function g 2 /(Aitr). 

(iii) The interaction vertex is ry°. 

For the simplicity of the notations, one can perform calculations using the unit system where 
v = 1, then restore v at the end results by dimensionality. To have analytic control, we shall 
work in the large-A limit and then extrapolate to iV = 2. Each fermion loop comes with 
a factor of N, so in the large- N limit one has to resum the fermion loop in the photon 
propagator (this is identical to the random phase approximation). The fermion loop is 
naively divergent; however, in any gauge-invariant regularization scheme (e.g., dimensional 
regularization) it is convergent. The resummed Coulomb propagator is (see the Appendix) 

Since D(q) ~ 1/N, even in the g 2 — > oo limit the interaction between fermions remains weak 
at large N, enabling a perturbative calculation. 



III. RENORMALIZATION GROUP AND THE STRONG-COUPLING FIXED 

POINT 

We shall now perform Wilson RG in our theory at the leading nontrivial order in 1/N. We 
assume that the theory has a cutoff Ao, and proceed to integrate out all momenta between 
Ai and Ao, where Ai is smaller than Ao by an exponential factor. The leading 1/N correction 
to the fermion kinetic term comes from the one-loop fermion self-energy graph, 

2 f d 3 q 7 (/-^)7° / g 2 N \q\ 2 . 



(2tt)3 (p - q y 




5 



where integration is over q in the momentum shell Ai < q < A . In RG, we are interested 
only in the contribution proportional to ln(Ao/Ai). 

Although the integral can be computed in closed analytic form for any A (see below), it 
is instructive to consider first the limit A = oo. In this limit, we find 

E(p) = S 7oPo + S17 • p, (11) 

where S and Ei are represented as integrals, 

y _ § f d*q q 2 -<? 8 f d*q qj 

N J (2tt) 3 {q 2 fl 2 \q^\ ' 1 N J (2vr) 3 (g 2 ) 3 / 2 |<f I ' 1 ' 

The integrands scale as g -3 ; therefore, both integrals contain the factor ln(Ao/Ai). These 
correspond to logarithmic renormalization of the kinetic terms ip^doip and 'tjrfdiij). The 
fact that these terms are renormalized differently means that the fermion velocity changes 
under the RG. 

However, it is easy to see that So and Ei contain an additional logarithmic divergence 
due to the singularity of the integrands in the limit \q\/qo — > 0. This singularity can be 
traced to the fact that the Coulomb interaction is unscreened in the finite frequency, zero 
wave number limit. However, in this limit the only effect of the gauge field is to phase rotate 
the fermion operator tp; therefore, the logarithmic singularity associated with the |^|/<7o ~~ * 
limit should disappear in any gauge-invariant quantity, for example, in the fermion velocity 
v . Indeed, the renormalization of v depends on Ei — So, which is free from the q = 
singularity: 

d 3 q 1 4 A 



Sl ~ S ° " N J J2^ WW 2 ~ 111 AT • (13) 
The RG equation for the velocity becomes 

= ~^ v(v) ' (14) 

which implies that in the limit of infinitely strong Coulomb coupling A —>■ oo, the velocity 
has a finite anomalous dimension 7^ = — 4/(7r 2 A^). The solution to the RG equation (1141) is 

v(p) = const x p - 4/(7r2jv) . (15) 

Since the velocity is the slope of the dispersion curve, one concludes that the fermion dis- 
persion relation has the form 

u) = const x p z , (16) 
with the dynamic critical exponent z being 

z = l + lv = l-J- + 0(N- 2 ). (17) 

7T z iV 



Note that the argument of Ref. [T7J that requires z = 1 does not apply to our case, since 
the fixed point here is at infinite Coulomb coupling g — > 00. 
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Some physical consequences of Eq. (fTTI) need to be mentioned. Since z < 1, the quasipar- 
ticle is stable, since its decay into two or more other quasiparticles is forbidden by energy and 
momentum conservation. The specific heat has a power-law behavior at small temperature 
T: 

C(T)~T a , a=--l. (18) 
At large N, a = 1 + 8/(7r 2 iV). The first logarithmic correction to the specific heat was 



computed in Ref. [18||; our formula ([TBI sums up all powers of log(T)/iV. 

The power-law behavior of the velocity ([TBI) tells us that at g — > oo the system is scale 
invariant with respect to the scaling transformation Q. In this regime, it is more convenient 
to define the dimensions of the operators with respect to the scale transformation ([3]) rather 
than the relativistic version (j3J). In this new scheme, 

N = -l, [t] = -z, [A ] = z, [g- 2 ] = l-z>0. (19) 

The last equation means that the bare kinetic term for A is a relevant perturbation at the 
strongly coupled fixed point. 

IV. FINITE COULOMB INTERACTION 

Let us briefly consider the case of finite A. The loop integral in Eq. (flQj) can be evaluated 
explicitly (see the Appendix), 

S(P) = ^[/o(A)po7° + A(A)p-7]ln^, (20) 

where 

A = ^ (2!) 

is the parameter measuring the importance of photon self-energy compared to its bare inverse 
propagator. The functions fo and f\ in Eq. (1201 are 

arccos A — 1 H — - , A < 1, 



A(A) = < _1_ 2A (22) 

VA " 1 lnfA + VA^T)-l + ^, A>1, 



A V / 2A 

2 - A 2 . vr x n 

arccos A — 2 + — , A < 1 , 



/o(A)=< , , A _ (23) 



= \n(\ + VX^t)-2 + j, A>1. 
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The asymptotics of the functions fo and /i at large and small A are: 

71 . A 2 



A (A) 



/o(A) 



4 



-A 



+ 0(A d 



TT 



ln(2A)-l + ^ + 0(A- 2 ) 
y + 0(A 3 ), 



7T 



ln(2A) - 2 + ^- + 0(A- 2 ) 
A 



A < 1, 
A > 1, 
A < 1, 
A > 1. 



(24) 



(25) 



These functions are related by /o(A) = — A /{(A) + /i(A) and are plotted in Fig. [TJ together 
with the difference fx — f$. The expressions above are consistent with those quoted in 
Ref. 
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FIG. 1: The functions /i(A) (solid line) and /o(A) (dashed line) and their difference /i(A) — /o(A) 
(dotted line). 



The running of the fermion velocity is governed by the RG equation 



P- 



dv{p) 
dp 



7T 



[/i(A)-/o(A)Mp)= 7 «(AMp). 



(26) 



Since /i(A) > f\X) for all A (see Fig. [TJ the fermion velocity increases monotonically as one 
decreases the momentum scale. Here, j v is the anomalous dimension for the velocity v. As 
expected, in the limit of infinite Coulomb coupling it approaches the constant value found 
above: 



lim 7„(A) 

A— >oo 



4 



TT 2 N 



(27) 



V. FINITE N AND GRAPHENE 



So far, we have discussed the limit ^> 1 where reliable calculations can be performed. 
Let us now discuss finite N, keeping in mind that in graphene N = 2. Unfortunately, there 
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is no reliable calculation tool that works for all values of N and the coupling constant g; 
therefore, our discussion will be mostly conjectural. 

The simplest possibility is that what we saw at large N remain qualitatively valid at all 
N: there are two fixed points, an infrared unstable one at g = oo and an infrared stable 
one at g = 0. At g = oo, the system exhibits invariance with respect to Eq. (jHJ), although 
the value of z for small iV cannot be computed in a reliable fashion. This possibility is 
illustrated in Fig. [2j If this is the case, then real graphene (with iV = 2) is always in the 
semimetal phase. 

00 , 



V 



semimetal 



N 



FIG. 2: The simplest phase diagram. The line g = oo is a line of infrared unstable fixed points. 
From any finite value of g, the system flows to an infrared stable fixed point at g = 0. 



Another possibility is that at sufficiently small N, the Coulomb interaction is strong 
enough to induce a spontaneous condensation of particle-hole pairs, creating an excitonic 
gap which makes the system insulating. (The alternative possibility of ferromagnetism was 



considered in Refs. [19|, |20j.) This possibility is depicted in Fig. [3j The insulator phase 
exists when N < A crit , but only for sufficiently strong coupling g > g c (N). When N = iV cr it, 
g c = oo, and for iV > A crit , the insulator phase no longer exists; the system is in the 
semimetal phase for all g. Relativistic massless (2+l)D QED is thought to develop a gap 
when number of fermion species is below some critical value [141 ]. 

If the phase diagram is as in Fig. [3j then two possibilities exist for real graphene. If 
A" crit < 2, or if N crit > 2 but the bare Coulomb coupling (in vacuum) g < g c (N = 2), then 
graphene is a semimetal. In contrast, if A cr ; t > 2 and in vacuum g > g c (N = 2), then 



freely suspended graphene [10| is an insulator. All available experimental data, on the other 
hand, are consistent with graphene on a Si02 substrate being a semimetal. Therefore, if in 
vacuum graphene is insulating, then it undergoes an insulator-semimetal phase transition 
as a function of the dielectric constant of the substrate. 
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insulator 



semimetal 



N 



N 



FIG. 3: A slightly more complicated phase diagram. For N > -/V cr it, the system is always in the 
semimetal phase. For N < N CT it, it can be in the semimetal phase [for g < g c (N)] or in the 
insulating phase [for g > g c (N)]. 



The authors of Refs. 



21 



22 



solved a gap equation with the screened Coulomb interaction 
and found N CT a ~ 2.55. If this is the case, then the system with N = 2 develops an excitonic 
gap at sufficiently large g. However, this result is probably not conclusive as the gap equation 



in Ref. 2l|, |22j is not systematic at small N, and also neglects the frequency dependence in 



the photon propagator. Ultimately, the phase diagram of graphene should (and could) be 
determined by a direct numerical simulation of the theory. 

In the rest of this section, we shall assume that the phase diagram is as in Fig. [2J or 
as in Fig. [3] with Nait < 2. The system with N = 2 and infinite Coulomb coupling should 
show scale invariance ([3]) with a critical exponent z which can be estimated by extrapolating 
Eq. (HZD, 

z « 0.8 for N = 2. (28) 

Note that the 1/N correction to z is only 20% at iV = 2, giving us hope that the 1/N 
expansion works reasonably well even for iV = 2. 

The question now is whether in real graphene A is large enough so that the system is closed 
to the scale-invariant regime. If we take the typical experimental value v = 10 6 m/s 0, 0] 
then we find A ~ 3.4, which is reasonably large compared to 1. Therefore, one concludes that 
graphene is sufficiently close to the scale-invariant regime. However, in most experiments 
the graphene single layer lies on a substrate, which reduces the effective Coulomb potential 
[Eq. (jBJ)]. If we take the substrate to be Si0 2 with e = 4.5, we find A ~ 1.25, which is only 
marginally larger than 1. 

The above quoted value A ~ 3.4 should be thought of as the initial condition for the RG 
equation, valid at the momentum scale of the order of the inverse lattice size. As we flow to 
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the infrared, the system deviates more and more from the strongly coupled fixed point due 
to the bare kinetic term which is a relevant perturbation. However as the initial condition is 
relatively close to the fixed point, it will take a large "RG time" to go to the weak coupling 
regime A < 1. The graphene thus remains close to the strong coupling limit for a large 
momentum range. The width of this range can be roughly estimated as 

exp — — In 3.4 ~ 10 3 , (29) 



where we have used the leading 1/N result for the dimension of the bare Coulomb kinetic 
term. 

An additional information for the effect of finite A can be obtained from the value of the 
anomalous dimension for the velocity as given by Eq. (|T4l) at A = 3.4: 

7,(3.4) « 0.7 7u (oo) « 0.15. (30) 

It is a 30% reduction of the anomalous dimension r y v . In Fig. [IJ we see that for the realistic 
A the different f\ — fo is already a rather flat function of A; an approximate scaling behavior 
can be expected. 

At the asymptotic infrared end of the RG flow is the trivial fixed point A = 0, near which 
the fermion velocity increases linearly with the logarithm of the momentum [as followed 
from Eq. (I14p and the small A asymptotics in Eqs. 



-P^ = A (31) 
1 dp 16tt' K ' 

which translates into an increase of 1.26 x 10 6 m/s per decade for graphene in vacuum. 
However, this regime is achieved only at extremely low momenta after the system has passed 
from the vicinity of the strongly coupled fixed point to the trivial one. 



VI. CONCLUSION 

In this paper, we consider a model of N massless four-component Dirac fermions interact- 
ing through an instantaneous Coulomb interaction. We show that in the limit of infinitely 
strong coupling, the system shows a scale-invariant behavior (J3]) which is characterized by 
a dynamic critical exponent z. We obtain the value of z at large N. 

We also discuss two possibilities for the phase diagram of the system at finite N . In one 
of the possibilities, a part of the phase diagram is occupied by the insulator phase. It is very 
interesting if in freely suspended graphene the Coulomb interaction is strong enough to open 
an excitonic gap. Provided that such a gap is never opened for iV = 2, we argue that the 
Coulomb interaction in vacuum is strong enough so that freely suspended graphene sheets 
can be considered as being in the vicinity of the strong-coupling fixed point for a large range 
of momentum. 
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How can our predictions be tested in experiment? Assume one could measure with high 
precision the quasiparticle dispersion curve in freely suspended graphene. If the strong- 
coupling large- N limit is any guide, then one should have a slight deviation from the linear 
dispersion law, perhaps approximately uo ~ p°' S5 [see Eq. (1301) ]. For graphene on a substrate, 
the deviation from the linear law is smaller. 

Finally, we hope that this work will motivate numerical simulations of the model (j3j), 
which should help clarify it phase structure. 
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APPENDIX A: CALCULATION OF THE BETA FUNCTION AT LARGE N 



1. Coulomb propagator 



To leading order in the large- N expansion, we must resum all fermion bubble graphs in 
the Coulomb propagator. We shall perform calculations in the unit system where v — 1 and 
restore v in final formulas when needed. 




FIG. 4: The photon polarization diagram. 



The one-loop fermion bubble diagram (Fig. H]) is 

^-^/IW^j^f?)- (Ai) 

where X — 7^m> M = 0, 1, 2. The manipulation of the Dirac algebra proceeds in a standard 
fashion. We use the formulas 

J = F- < A2 > 

Tr (7V7V) = 4(« pCT - 5^5™ + 5^5 up ) (A3) 
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to find 

TV ^ a at 2 f d3k 2k o( k ° + Qo) -k-(k + q) 

We now use Feynman parametrization 

1 f 1 1 

AB = J dX [xA + (1 - x)B] 2 (A5) 

to rewrite 

n (5 ) = iNf f dx {** ph^Lz^L^ . (A6) 

Changing the integration varibale k — > A; — xg, the denominator in the integrand becomes an 
even function of q, and one can throw away terms odd in k in the numerator. Furthermore 
due to spherical symmetry we can replace in the numerator k% —>■ k 2 /3. We find 

n(9) = 4^ \\ -^-^:^:^- x)f . ( A7) 







(2tt) 3 [k 2 + x(l - x)q 

We perform integration over k using standard formulas of dimensional regularization 



d d k 11 r(n-|) 1 



(2yr) d (fc 2 + A) n (47r) rf / 2 r(n) A"" d / 2 
d d k k 2 1 dr(n-f-l) 



(2vr) rf (k 2 + A) n (4n) d / 2 2 T(n) &n-d/2-i 



(A8) 
(A9) 



n(g) = ^^fd* [x(i - = , (Aio) 



to obtain 

irr^ _ * _ .... . . . , , 

n J 8 \/qJTv 2 \q] 

where in the last expression we have restored v. The resummed Coulomb propagator is 

DU,)=lm + 9 -T^=\ (All) 



2. Correction to fermion propagator 

We now compute the 1/N correction to the fermion self-energy (Fig. [5]): 

= ~ g jj^f (p-q) 2 ° {q) - (A12) 

The integral is naively linearly divergent at large q, but to leading order in q the integrand 
is an odd function of q. Therefore, the integral is only logarithmically divergent and can be 
evaluated by expanding in p <C q. One finds 

E(p) = Z ol °p + Z 1 Yp ll (A13) 
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FIG. 5: The fermion self-energy diagram. 



where 



d 3 q ql - \q\ 



2 



Zo = 9 2 77T^ JL -T^D(q), (A14) 



(2tt) 3 q 

z ^s 2 Jwr4 Diq) - (A15) 

Introducing spherical coordinates: go — qcos8, |g| = qsind, one can then write 

4 A J . n ,„ cos 2 6- sin 2 6 fdq , k . 

- sin6 d6 - , / — , (A16) 



ir 2 N 2 J sin 0(1 + A sin i 



o 



4 A f . n in cos 2 9 fdq , . 1 nS 

Z x = — sm6de ; r/— , (Al7) 

1 n 2 N 2 ./ sin0(l + Asin0) J q ' V ; 



o 

,2 ; 



where A = g N/16. The integral over tig in a spherical shell A x < g < A yields ln(Ao/Ai] 
The angular integral can be taken explicitly and results in Eqs. fl20|) . fl22l) . and (1231) . 
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